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(4)  Introduction 

The  ability  of  cancer  cells  to  proliferate  inapproriately  is  due  the  loss  of  tumor 
suppressor  genes  and  the  gain  of  function  of  oncogenes.  Two  classes  of  tumor  suppressor 
genes  exist,  mutators  and  growth  regulators.  Mutators  are  genes  that  when  mutant  cause  an 
increase  in  the  rate  of  genomic  instability  and  hence  the  rapid  accumulation  of  mutations  of 
all  classes.  Genes  of  the  mutator  class  include  mis-match  repair  genes  involved  in  HPCC, 
DNA  repair  genes  such  as  X.P.,  cell  cycle  checkpoint  genes  such  as  ATM  and  p53. 
Recently  two  breast  cancer  genes  BRCA1  and  BRCA2  have  been  implicated  in  DNA  repair 
and  fall  into  this  class.  The  second  class  of  tumor  suppressors,  growth  regulators,  are  those 
directly  involved  in  regulating  cell  proliferation  or  the  ability  of  tumor  cells  to  survive  and 
metastasize.  Oncogenes  such  as  myc  or  ras  are  dominant  and  act  in  opposition  to  tumor 
suppressors.  Many  oncogenes  and  tumor  suppressors  regulate  the  cell  cycle  and  affect  celll 
cycle  regulation.  In  order  to  understand  this  regulation,  it  is  important  to  have  knowledge  of 
genes  under  cell  cycle  transcriptional  regulation.  Mitosis  represents  one  of  the  most 
fundamental  activities  of  eukaryotic  cells,  playing  an  essential  role  in  events  ranging  from 
gametogenesis  and  wound-healing  to  multicellular  development1'-3.  Proper  regulation  of 
gene  activity  during  the  cell  cycle  is  therefore  likely  to  govern  critical  aspects  of  diverse 
biological  processes.  Indeed,  disruption  of  cell  cycle-dependent  mRNA  regulation  results 
in  profound  phenotypic  consequences.  For  example,  defects  in  the  retinoblastoma  (Rb)  and 
E2F-1  proteins,  which  control  transcription  at  the  G,/S  transition,  precipitate  developmental 
defects,  aberrant  cell  morphologies,  and  uncontrolled  proliferation4'6.  The  pleiotropic 
effects  of  these  mutations  suggest  that  a  considerable  proportion  of  all  genes  experience 
differential  regulation  during  the  cell  cycle.  However,  although  these  transcripts  have  been 
comprehensively  catalogued  in  yeast,  few  targets  of  cell  cycle-dependent  transcriptional 
machinery  have  been  identified  in  higher  organisms  A  .  A  global  perspective  on  genetic 
regulation  during  eukaryotic  mitosis  might  therefore  significantly  expand  our  understanding 
of  both  normal  cell  division  and  the  abnormal  phenotypes  observed  in  certain  pathological 
processes. 

In  addition,,  the  transcriptional  induction  of  genes  in  response  to  DNA  damage  is  a 
critical  component  of  the  cellular  response  to  DNA  damage  and  the  prevention  of 
mutagenesis.  The  identification  of  transcriptional  target  genes  is  critical  to  understanding 
how  these  regulatory  pathways  function.  In  this  grant  we  propose  the  identification  of 
target  genes  such  as  these  to  better  our  understanding  of  these  processes  and  to  provide 
tools  for  future  analysis. 

(5)  Body 

As  noted  in  our  previous  reports,  during  the  process  of  devising  retroviral  vectors  to 
identify  genes  under  the  control  of  promoters  regulated  by  DNA  damage  and  other  stimuli, 
we  began  a  collaboration  with  several  groups  to  identify  genes  induced  by  these  stimuli 
using  transcriptional  arrays.  Our  retrovial  work  was  never  able  to  identify  more  than  just  a 
few  genes  that  were  reproducibly  induced  by  DNA  damage.  Thus  due  to  the  massive 
amount  of  work  involved  in  identifying  the  gene  for  each  insertion,  we  opted  to  perform 
transcriptional  profiling  using  microarrays.  This  report  details  our  success  in  this  area. 

Transcriptional  Profiling  of  Cell  Division  and  DNA  Damage 

In  this  study,  we  measured  cell  cycle  regulated  transcripts  in  human  cells  on  a 
genome-wide  scale  using  high-density  oligonucleotide  arrays.  Primary  fibroblasts  prepared 
from  human  foreskin  were  grown  to  approximately  30%  confluence  and  synchronized  in 
late  G,  using  a  double  thymidine  block  protocol9.  Cultures  were  released  from  arrest,  and 
cells  were  collected  eveiy  two  hours  for  twenty-four  hours,  covering  nearly  two  complete 
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cell  cycles.  Messenger  RNA  was  isolated,  labelled,  and  hybridized  to  arrays  containing 
probes  for  approximately  40,000  human  genes  and  non-overlapping  ESTs^.  The  entire 
experiment  was  performed  in  duplicate  under  identical  conditions  for  a  subset  of  these 
genes  for  an  independent  confirmation  of  results. 

The  two  cell  cycle  data  sets  were  concatenated  and  analyzed  using  both  supervised 
and  unsupervised  clustering  of  expression  patterns  (See  Methods).  Five  groups  of 
transcripts  were  identified  that  displayed  a  periodicity  consistent  with  the  length  of  the  cell 
cycle  in  both  data  sets,  with  peaks  corresponding  to  early  G,,  late  G,,  S,  G?,  and  M  phase. 
Using  this  approach,  1,332  transcripts,  including  895  previously  uncharacterized  ESTs,  were 
assigned  to  cell  cycle-regulated  expression  clusters.  These  assignments  were  made  only  for 
transcripts  that  clearly  oscillated  in  obvious  patterns,  and  the  conservative  nature  of  this 
analysis  makes  it  possible  that  some  cell  cycle  regulated  transcripts  were  not  selected. 
Consistent  with  previous  studies,  multiple  cyclin  family  members  were  identified,  as  well  as 
Cdc2,  Cdc8,  Cdc25  and  numerous  genes  involved  in  DNA  replication  and  chromatin 
structure^.  This  information  can  be  viewed  and  downloaded  at  the  Internet  address 
www.mag.com/expressionanalvsis. 

We  also  sought  to  better  understand  the  relationship  between  transcriptional 
changes  that  normally  occur  during  mitosis  and  those  observed  during  events  that  affect  the 
progression  of  cell  division.  For  the  purposes  of  this  assessment,  transcript  levels  were 
quantitated  in  fibroblasts  after  exposure  to  ultraviolet  light  and  methyl  methane  sulfonate 
(MMS),  both  of  which  precipitate  DNA  damage  and  consequent  cell  cycle  arrest  at  multiple 
checkpoints. 

All  experiments  were  performed  in  duplicate  and  compared  with  untreated  cells  to 
identify  consistent  differences.  Stringent  requirements  were  also  used  to  identify  these 
differences.  To  qualify,  changes  had  to  be  called  in  at  least  two  independent  experiments 
and  more  than  a  two-fold  difference  was  required.  Of  the  more  than  6,000  genes  assayed, 
64  were  consistently  induced  and  24  repressed  after  UV  exposure,  and  73  were  induced  and 
18  repressed  after  MMS  exposure.  These  data  are  also  available  at  the  Internet  site  given 
above. 

Systematic  Functional  Analysis  of  Expression  Clusters 

Previous  large-scale  expression  studies  have  relied  heavily  on  subjective  evaluations 
of  the  functions  of  differentially  expressed  genes^*  11,12  However,  the  rapidly  increasing 
number  of  these  data  sets  has  created  an  urgent  need  for  more  standardized  methods  for 
assessing  the  coordinate  regulation  of  biological  pathways.  Overrepresentation  of  genes 
from  a  common  functional  class  within  an  expression  cluster  may  reveal  global  pathway 
activation^.  We  implemented  a  statistical  analysis  based  on  this  hypothesis.  More  than 
2400  characterized  human  genes  represented  on  the  arrays  were  classified  into  160 
functional  categories  (Methods).  Using  the  binomial  distribution  function,  P  values  were 
calculated  based  on  the  frequency  of  occurrence  of  genes  with  common  biological  function 
in  each  expression  cluster.  Statistically  significant  enrichments  for  specific  functions  were 
detected  in  four  of  the  five  cycle-regulated  transcription  clusters,  with  P  values  ranging  from 
1.4  x  10‘4  to  1.9  x  10". 

Transcriptional  Regulation  of  Central  Mitotic  Processes 

The  most  basic  test  of  the  systematic  analysis  applied  in  this  study  is  the  detection 
of  biological  processes  known  to  play  a  central  role  in  cell  division.  As  anticipated,  we 
found  coordinate  upregulation  of  transcripts  involved  in  cell  cycle  control,  DNA  replication, 
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chromosome  segregation,  and  cytokinesis,  providing  an  important  validation  of  this 
approach. 

The  comprehensive  identification  of  specific  genes  that  experience  cell  cycle- 
dependent  transcriptional  control  in  these  known  pathways  may  provide  insight  into  how 
regulatory  disruption  may  exert  biological  effects.  We  first  examined  transcripts  for  genes 
that  govern  cell  cycle  progression  (P  =  4.3  x  1012)  Cyclins  such  as  A,  B,  and  F  family 
members  were  coordinate^  induced  at  the  G2-M  transition,  at  which  point  they  regulate  the 
entry  into  mitosis.  However,  some  genes  that  inhibit  the  G,-S  transition  (such  as  pl6 INK4 
and  Mytl)  were  induced  at  this  time  period,  while  others  (such  as  pllKipl)  were  maximally 
repressed.  The  most  striking  functional  coordination  was  observed  for  genes  involved  in 
proteolytic  control  of  cell  cycle  progression,  including  Cdc5,  Cdc20,  Cksl,  Cdc34,  E2A,  and 
the  cyclin-specific  ubiquitin  ligase  UbclO,  all  of  which  are  induced  during  M  phase. 
Notably,  a  second  set  of  ubiquitin-associated  transcripts  was  highly  expressed  during  the 
GrS  transition,  including  the  Nedd4-like  ubiquitin  protein  ligase  Wwp2,  ubiquitin-specific 
protease  Usp7,  and  cullin  Cul4A. 

Consistent  with  findings  in  yeast,  numerous  DNA  replication  genes  were  induced  in 
late  G,  and  S  phase,  including  some  genes  known  to  be  positively  regulated  by  E2F 
transcription  factors,  such  as  DNA  polymerase  alpha,  DNA  primase,  components  of 
replication  factor  A  and  C,  MCM  genes,  histones,  and  dihydrofolate  reductase  ( P  =  4.4  x 
10'5)  (Table  2).  Interestingly,  we  also  observed  unexpected  overrepresentation  of  genes 
related  to  actomyosin  function  in  the  S  phase  cluster  (P  =  6.7  x  10  4).  Many  of  these 
transcripts  encode  proteins  directly  implicated  in  cytokinesis,  including  several  smooth 
muscle  myosins,  alpha-actin,  the  small  GTPase  Rho,  gelsolin  Sm22,  calponin  Cnnl,  and 
Tau.  Induction  during  S  phase  was  also  observed  for  numerous  effectors  of  actin-based 
cytoskeletal  remodeling  not  previously  linked  to  cell  division;  for  example,  the  regulatory 
genes  Arf6,  Cdk5,  Tiaml,  and  Hefl.  Other  transcripts  found  in  this  cluster  encode  gelsolin¬ 
like  actin  polymerization  domains  (such  as  Flil  and  Capg ),  but  have  not  been  assigned  a 
more  specific  biological  role. 

Genes  that  regulate  chromosome  segregation  were  coordinately  upregulated  during 
M  phase,  including  the  spindle  checkpoint  sensor  Mpsl,  the  centrosomal  kinase  Nek2,  and 
stathmin,  a  putative  control  point  for  microtubule  assembly  (P  =  8.4  x  lO10)1'*'1”.  Nearly 
global  upregulation  during  M  phase  was  also  observed  for  microtubule-dependent  motor 
proteins  that  physically  effect  chromosome  segregation,  such  as  mitosin  and  virtually  every 
kinesin  known  to  interact  directly  with  the  centromere1'’111.  Qear  cen  CyC]e  regulation  was 
not  observed  for  dynein  genes,  suggesting  restriction  of  this  transcriptional  control  to  a 
subset  of  motor  proteins. 

Cell  Cycle  Control  of  Extracellular  Matrix  Activity  and  Apoptosis 

These  analyses  reveal  that  an  unexpectedly  high  proportion  of  the  genes  induced 
during  G2  phase  function  in  the  extracellular  matrix  (P  =  1.9  x  10  ").  Genes  found  in  this 
cluster  encode  cell-cell  adhesion  molecules  such  as  collagen  VI  alpha- 1,  integrin  alpha-6, 
desmoplakin,  plakoglobin,  and  laminin.  Other  matrix-associated  genes  found  in  both  the  G2 
and  M  clusters  positively  regulate  cell  movement  through  the  extracellular  matrix;  for 
example,  platelet-derived  growth  factor  receptor  PdgfA,  vasodilator-stimulated 
phosphoprotein  Vasp,  hyaluron  synthase  Has2,  hyaluron  receptor  Rhamm,  tissue 
plasminogen  activator  Plat,  and  urokinase  plasminogen  activator  receptor1 9-22, 

Marked  enrichment  was  also  detected  in  the  late  S  phase  cluster  for  Bel-related 
transcripts  (P  =  2.1  x  10'4),  which  have  been  extensively  characterized  as  inhibitors  of 
apoptotic  proteases23.  These  genes  encode  Bcl-7A,  Bcl-xL,  Bak,  Mcl-1,  the  Bcl-2 
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interacting  proteins  Bbc3  and  Bad,  and  the  Bcl-xl  interacting  protein  Bap31.  Interestingly, 
no  other  group  of  apoptotic  inhibitors  were  overrepresented  in  any  of  the  cell  cycle 
expression  clusters.  It  was  also  observed  that  none  of  the  Fas  or  TNF-related  genes  that 
positively  regulate  apoptosis  were  found  to  show  a  significant  degree  of  cell  cycle  regulation 
in  these  data. 

DNA  Damage  Response  Involves  S  phase  Transcripts 

It  is  possible  that  cell  cycle  regulated  transcripts  comprise  organized  programs  of 
biological  activities  that  are  required  outside  of  cell  division.  DNA  damage  in 
Saccharomyces  cerevisiae,  for  example,  is  known  to  induce  a  significant  proportion  of  the 
transcripts  induced  during  late  G,  phase  (M.  Campbell  and  S.  Tavazoie,  unpublished  data). 
Such  overlapping  specificities  may  identify  cellular  responses  that  involve  cell  cycle- 
dependent  transcriptional  machinery. 

In  the  6,800  transcripts  surveyed  in  this  study,  42  were  reproducibly  induced  in 
response  to  both  UV  light  and  MMS  in  human  fibroblasts.  This  group  included  many 
known  damage-induced  transcripts  such  as  c-FOS,  JUNB,  GADD45,  and  IEX1.  Strikingly, 
these  genes  were  found  to  share  considerable  target  specificity  with  transcripts  in  the  S 
phase  cluster  (P  =  2.3  x  10'8)  but  not  with  any  other  group  of  cell  cycle  regulated  genes. 
FACS  analyses  indicate  that  these  cells  were  not  arrested  in  S  phase  when  samples  were 
collected,  excluding  the  possibility  that  such  an  effect  was  responsible  for  the  observed 
induction  of  these  transcripts. 

Yeast  genes  induced  by  DNA  damage 

Since  the  goal  of  this  project  was  to  identify  DNA  damage  inducible  genes  from 
eukaryotes,  we  thought  it  would  be  useful  to  also  know  the  spectrum  of  DNA  damage- 
inducible  genes  in  a  model  eukaryote  S.  cerevisiae.  We  have  established  a  collaboration  with 
Pat  Brown's  lab  to  establish  the  complete  list  of  DNA  damage  inducible  genes  from  yeast 
using  chip  hybridization  technology9.  We  have  performed  a  number  of  analyses  and  have 
identified  approximately  300  transcripts  that  are  significantly  increased  in  response  to  DNA 
damage.  Many  of  these  have  obvious  human  homologs  that  we  can  test  for  DNA  damage- 
inducibility.  Interestingly  we  found  that  they  fall  into  several  categories  based  on  time  of 
induction  during  a  time  course  and  into  groups  of  genes  based  on  function.  These  genes 
are  shown  grouped  into  functional  units  in  Tables  1  and  2  below.  In  addition  to  genes 
involved  in  DNA  replication,  nucleotide  synthesis  and  DNA  repair  we  were  surprised  to  see 
a  number  of  genes  involved  in  electron  transport,  carbohydrate  metabolism,  Amino  acid 
biosynthesis,  protein  folding  and  degradation.  This  data  has  now  been  submitted  for 
publication. 
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Table  1 


Types  of  yeast  genes  induced  by  DNA  damage  by  MMS  through  a  2  hr  time  course. 

Total  number  of 

number  of  genes 

number  of  genes 

number  of  genes 

genes  induced  over 

peaking  at  15-30  min 

peaking  at  45-60  min 

peaking  at  90-120 

2-fold 

min 

Amino  Acid  Biosynthesis 

22 

14 

4 

4 

Protein  Degradation 
39 

6 

13 

20 

Protein  Folding 

7 

6 

1 

0 

Transport 

20 

7 

4 

9 

Nucleotide  synthesis 
15 

0 

1 

14 

Stress-induced  genes 
20 

3 

2 

15 

Electron  transport 

8 

0 

2 

6 

Carbohydrate  metabolism 

38 

45 

0 

7 

DNA  replication  and  repair 

13 

0 

6 

7 

Wild  type  cells  were  treated  with  0.03%MMS  (80%  viability) 
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Table  2.  Names  of  4  classes  of  yeast  genes  induced  by  DNA  damage 


DNA  replication  and/or  repair 

DUN  1  Checkpoint  protein  kinase 
RNR1*  Ribonucleotide  reductase 
RNR2*  Ribonucleotide  reductase 
RNR4*  Ribonucleotide  reductase 
TRX2* 

MAGI  Methyl  transferase 
RFA1  DNA  Replication 
RFA2  DNA  Replication 
POL30  DNA  Replication 
PHR1  Photo-repair 
RAD51  Recombination 
TRR1 

YHR029C  *  (similar  to  thymidylate 
synthase) 


Stress-induced  genes 
SOD1  Superoxide  dismutase 
HYR1 

YAP1  Transcription  factor 

YHB1 

TIR1 

XBP1  Transcription  factor 
TIR2 

HSP42  Heat  shock  protein 

HSP104  Heat  shock  protein 

HSP82  Heat  shock  protein 

HSP12  Heat  shock  protein 

HSP26  Heat  shock  protein 

SSE2  Heat  shock  protein 

STI1  Heat  shock  protein 

SIS1  Chaperone 

YGP1 

US  VI 

DDR48 

GRE2 

GRE3 


Nucleotide  Biosynthesis 

THI1 

THI5 

THI12 

TRR1 

TRR2 

UGP1 

YNK1 

ADE17 

STF2 

YBR284W  (similar  to  AMP  deaminase) 
RNR1* 

RNR2* 

RNR4* 

TRX2* 

YHR029C  *(similar  to  thymidylate 
synthase) 

Genes  involved  in  protein  degradation 

RPN6 

RPT4 

RPT6 

RPN7 

RPT2 

RPN12 

RPT1 

RPN3 

RPT3 

RPT5 

PRE5 

PRE9 

PRE10 

PUP3 

PRE3 

PRE2 

PRE8 

PUP1 

PRE1 

UMP1 

UBA1 

UBC5 

UBC8 

UBC13 

UBI4 

UFD1 

YPS6 

YPS6 

APE2 

PRB1 

PEP4 
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PRD1 

LAP4 

CDC48 

UBP6 

YAP180 

AUT7 

APG1 

APG7 


Discussion 

The  expression-based  analyses  used  in  this  study  revealed  both  known  and 
unexpected  results.  DNA  replication  genes  were  strongly  upregulated  during  late  Gp  as 
were  numerous  cell  cycle  control  and  chromosome  segregation  genes  during  M  phase 
(Table  1).  However,  the  observed  global  upregulation  of  genes  related  to  actomyosin 
function  and  programmed  cell  death  during  S  phase,  and  induction  of  extracellular  matrix- 
associated  genes  during  G2,  indicates  that  a  much  broader  set  of  biological  pathways  is 
transcriptionally  activated  during  cell  division  than  was  previously  known.  Notably, 
whereas  eight  biological  functional  classes  were  found  to  be  strongly  overrepresented  in  cell 
cycle  regulated  expression  clusters,  none  were  significantly  underrepresented.  One 
interpretation  of  this  observation  is  that  acquisition  of  cell  cycle-specific  regulatory  elements 
occurs  stochastically  in  the  genome,  but  is  retained  by  genes  for  which  mitotic  regulation 
confers  a  biological  advantage. 

These  data  offer  new  insights  into  established  mechanisms  of  cell  cycle  regulation. 
For  example,  the  coordinate  induction  of  ubiquitin-associated  cell  cycle  control  genes 
during  M  phase  indicates  a  requirement  for  restricting  the  time  of  activity  of  these  genes, 
and  may  reveal  how  aberrant  transcription  directly  interferes  with  proteolytic  cell  cycle 
regulation.  Interestingly,  these  results  are  in  direct  contrast  to  observations  made  in 
Saccharomyces  cerevisiae,  where  transcriptional  control  of  proteolytic  regulators  is 
observed  only  rarely^.  Recent  evidence  has  implicated  cullin  family  members  in  the 
proteolytic  degradation  of  G,-S  regulators,  supporting  the  possibility  that  the  ubiquitin- 
associated  genes  induced  at  this  time  may  participate  in  checkpoint  regulation^. 

These  data  also  implicate  previously  characterized  genes  in  specific  mitotic  activities. 
For  example,  the  nucleotide  excision  repair  genes  Rad2  and  Erccl  are  induced  during  S 
phase,  suggesting  that  they  may  function  in  the  correction  of  misincorporation  in  newly 
synthesized  DNA.  We  also  observed  that  the  DNA  damage-responsive  gene  Gadd45, 
increased  levels  of  which  have  been  shown  to  block  passage  through  a  G2-M  checkpoint,  is 
transcriptionally  repressed  during  this  period^.  Naturally  low  levels  of  Gadd45  transcript 
during  G2-M  may  permit  cell  cycle  progression  in  the  absence  of  DNA  damage-dependent 
induction,  providing  a  sensitive  means  of  checkpoint  regulation^. 

Examination  of  the  actomyosin-related  transcripts  observed  to  be  induced  during  S 
phase  reveals  that  many  of  these  genes  govern  assembly  of  the  contractile  ring.  Rho- 
mediated  reorganization  of  actin  stress  fibers  is  essential  for  cytokinesis,  and  its 
deregulation  has  been  implicated  in  tumor-specific  abnormalities  in  cell  separation^. 
Calponin  contains  protein  domains  that  are  required  in  yeast  for  actin  binding  during 
cytokinesis,  while  the  bacterial  homolog  of  Tau  has  been  shown  to  maintain  integrity  of  the 
contractile  ring^S, 29  Coordinate  transcriptional  upregulation  of  this  set  of  genes  occurs 
unexpectedly  early  in  the  cell  cycle,  well  before  the  onset  of  mitosis. 
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A  number  of  the  actin-related  transcripts  induced  during  S  phase  are  not  directly 
implicated  in  cytokinesis.  Their  gene  products  may  be  responsible  for  more  global 
cytoskeletal  changes  that  accompany  mitosis.  For  example,  Arf6  directs  reorganization  of 
cortical  actin  at  both  endosomes  and  the  plasma  membrane  and  has  recently  been  implicated 
in  cell  migratory  activities30>31.  Both  Cdk5  and  Tiaml  modulate  Rac  activity  and  are 
required  for  the  proper  development  of  actin-based  outgrowths  such  as  neurites^2.  Other 
genes  such  as  Flil  and  Capg  encode  gelsolin-like  actin  polymerization  domains  and  are 
hypothesized  to  function  in  cellular  motility,  but  have  not  been  assigned  a  specific  role  in 
cell  division^.  These  results  suggest  that  diverse  cytoskeletal  remodeling  pathways  are 
activated  to  facilitate  the  gross  morphological  changes  required  during  cell  division. 

The  observed  cell  cycle  regulation  of  these  numerous  cytoskeleton-related  signaling 
targets  offers  an  opportunity  to  examine  the  temporal  organization  of  their  transcriptional 
behavior.  The  induction  in  S  phase  of  Rho  with  regulatory  and  downstream  effector  genes 
such  as  Rho  GDP  dissociation  inhibitor,  actin,  and  Tau ,  and  the  simultaneous  expression  of 
Rac  with  genes  such  as  Atf6,  Cdk5,  and  Tiaml  suggest  that  genes  in  the  same  signaling 
pathway  are  generally  co-regulated.  We  examined  cell  cycle  regulated  genes  functionally 
related  to  another  small  Ras-like  GTPase,  Cdc42,  which  is  known  to  participate  in  signaling 
cascades  for  actin  and  filamin  reorganization.  Specifically,  Ack  kinase  and  components  of 
the  Arp2/3  actin  polymerization  complex  have  been  linked  to  signaling  pathways  involving 
Cdc42,  but  not  Rho  or  Rac^4,35.  Interestingly,  we  found  that  transcripts  encoding  all  of 
these  proteins  were  co-induced  during  G2  phase.  Such  temporal  compartmentalization 
could  increase  discrimination  between  closely  related  signaling  pathways.  It  is  also  possible 
that  phase  of  induction  separates  cytoskeletal  remodeling  programs  that  must  be  activated 
sequentially  during  mitosis. 

The  systematic  analysis  applied  here  is  perhaps  most  valuable  in  detecting 
unexpected  coordinate  regulation  of  functional  pathways,  which  is  likely  to  prove  difficult  to 
assess  based  solely  on  subjective  surveys  of  genome-scale  experimental  data.  We  identified 
an  overrepresentation  of  transcripts  with  extracellular  matrix-associated  functions  in  the  G2 
cluster.  At  least  some  of  these  gene  products,  such  as  plakoglobin,  desmoplakin,  integrin 
alpha-6,  and  laminin,  represent  components  of  structures  such  as  desmosomes  and 
hemidesmosomes,  which  facilitate  intercellular  adhesion  and  communication^. 
Transcriptional  upregulation  of  these  genes  likely  prepares  dividing  cells  for  re¬ 
establishment  of  contact  and  communication  with  the  extracellular  environment.  However,  a 
far  greater  number  of  matrix-associated  transcripts  in  the  G2  and  M  clusters  are  known  to 
positively  regulate  cellular  motility,  in  some  cases  through  signaling  to  the  actin 
cytoskeleton.  This  group  of  genes  includes  HAS2,  RHAMM,  VASP,  PLAT,  and  urokinase 
plasminogen  activator.  It  is  possible  that  these  genes  play  a  major  role  in  the  separation  and 
migration  of  daughter  cells  following  mitosis.  Interestingly,  these  gene  products  have  been 
characterized  most  extensively  in  terms  of  their  ability  to  increase  the  motility  of  tumor  cells. 
Overexpression  of  either  Has2  or  Rhamm  is  sufficient  to  promote  development  of  invasive 
cellular  phenotypes20>22  Therefore,  restriction  of  the  activity  of  these  genes  to  defined  cell 
cycle  intervals  may  represent  a  general  mechanism  for  enabling  activites  required  for 
proliferation  while  limiting  tumorigenic  cellular  behavior. 

The  global  repression  of  Bcl-2  related  antiapoptotic  genes  throughout  G2  and  M  is 
consistent  with  findings  that  deactivation  of  this  particular  class  of  genes  is  partially 
responsible  for  an  increased  susceptibility  to  programmed  cell  death  during  this  period^ ' . 
The  induction  of  antiapoptotic  genes  during  S  phase  may  also  serve  to  counteract  the 
proapoptotic  effects  of  the  E2F-1  transcription  factor,  which  positively  regulates  the  G,/S 
transition^.  Recent  evidence  has  supported  an  unexpectedly  broad  role  for  Bel-related 
genes  in  cell  cycle  progression,  possibly  through  direct  regulation  of  the  E2F  class  of 
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transcription  factors  or  other  activities  downstream  of  cyclin-dependent  kinase  2  (Cdk2)39- 
41.  One  candidate  effector  for  cell  cycle  regulation  of  5c/-related  genes  is  the  Stat5 
transcription  factor,  which  is  known  to  positively  regulate  transcription  ffomthe  Bcl-xL 
promoter  and  whose  expression  we  also  observed  to  be  restricted  to  S  phase42.  Further 
experiments  will  be  required  to  evaluate  the  full  functional  significance  of  these 
observations. 

The  significant  overlap  between  DNA-damage  induced  and  S  phase  transcripts 
indicates  that  common  subsets  of  genes  may  be  required  for  response  to  apparently 
unrelated  cellular  conditions.  This  phenomenon  has  been  observed  in  S.  cerevisiae,  where 
widely  divergent  biological  processes  have  been  linked  by  their  dependence  on  the  proper 
function  of  a  common  set  of  genes43.  One  explanation  might  be  that  numerous  genes  that 
function  directly  in  DNA  damage  repair  are  required  for  the  proper  execution  of 
chromosome  replication.  However,  none  of  the  transcripts  that  appear  in  both  groups 
participate  directly  in  DNA  synthesis.  Instead,  these  genes  function  predominantly  in  the 
stress  response,  either  through  intracellular  signaling  (the  transcription  factors  Jun-B,  Atf3, 
Atf4,  Egrl,  and  Egr2  and  the  dual  specificity  phosphatases  Duspl  and  Dusp5 )  or 
extracellular  communication  (ligands  including  IL-6,  IL-11,  and  EGF- like  growth  factor). 
Therefore,  it  appears  likely  that  the  general  regulatory  response  to  cellular  stress  is  governed 
partially  by  transcriptional  mechanisms  that  also  coordinate  gene  activity  during  the  cell 
cycle.  It  remains  to  be  determined  what  other  transcriptional  responses  display  significant 
similarities  with  the  expression  clusters  identified  from  these  data. 

A  number  of  genes  that  are  clearly  cell  cycle-regulated,  including  Edn-1,  Cox-2, 
Tfpi-2,  Plat,  Pai2,  and  Plaur,  were  previously  reported  to  be  transcriptionally  activated 
effectors  of  the  serum  response^.  Because  one  effect  of  serum  deprivation  is  mitotic 
synchronization,  the  previously  observed  transcriptional  results  are  unlikely  be  specific  to 
the  serum  response.  Further  analysis  and  comparison  with  the  previous  results  will  allow 
identification  of  serum-dependent  transcriptional  effects  that  are  not  simply  the  result  of  cell 
cycle  arrest.  These  data  should  provide  a  helpful  general  reference  point  for  future  human 
profiling  studies  that  involve  aspects  of  cell  cycle  progression. 

More  than  800  transcripts  from  previously  uncharacterized  EST  sequences  were 
also  found  to  be  clearly  cell  cycle  regulated.  Because  many  human  open-reading  frames 
have  no  significant  similarity  to  characterized  genes,  prediction  of  biological  function  based 
entirely  on  sequence  homology  is  limited.  Identification  of  transcriptional  behavior  under 
different  conditions  provides  an  additional  means  for  assessing  the  function  of  these  genes. 
Approximately  60%  of  the  previously  characterized  genes  identified  in  this  study  as  cell 
cycle  regulated  were  known  to  play  a  mitotic  phase-specific  role,  compared  to  5%  of  all 
genes  represented  on  the  arrays.  Therefore,  observation  of  cell  cycle  regulation  for  a 
transcript  indicates  a  more  than  10-fold  increase  in  the  probability  of  mitotic  function. 

The  cell  cycle  cluster  to  which  an  uncharacterized  EST  is  assigned  provides 
additional  clues  to  its  biological  role.  For  example,  based  on  the  calculations  made  about 
biological  function  in  this  paper,  genes  upregulated  during  M  phase  have  a  more  than  20- 
fold  higher  likelihood  than  random  genes  of  participating  in  cell  to  cell  communication. 
Furthermore,  assessing  these  transcriptional  results  in  the  context  of  sequence  homology 
may  provide  additional  clues  to  cellular  rolell.  For  example,  one  uncharacterized  EST 
induced  during  G2/M  shows  similarity  to  katanin.  an  ATPase  whose  oligomerization  leads 
to  microtubule  severing  during  cell  division44,45.  Another  EST,  induced  in  S  phase, 
displays  high  amino  acid  identity  with  the  RTEF1  transcription  factor,  which  preferentially 
induces  muscle-specific  genes4o,47  This  gene  may  therefore  be  partially  responsible  for 
the  global  cell  cycle  regulation  of  actin  and  myosin  transcripts.  The  increasing  number  of 
genome-scale  transcriptional  data  sets  will  soon  generate  multiple  lines  of  evidence  about 
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the  cellular  role  of  many  newly  sequenced  or  previously  uncharacterized  genes.  These  data 
should  provide  opportunities  for  more  systematic  comparisons  between  expression  and 
DNA  sequence  data  for  the  purposes  of  predicting  gene  function. 

Although  the  perspective  provided  by  this  study  is  restricted  to  the  transcriptional 
level,  the  detection  of  functions  ranging  from  DNA  replication  to  cytoskeletal  reorganization 
indicate  that  these  data  compose  a  representative  survey  of  mitotic  activities  in  higher 
eukaryotes.  The  global  regulation  of  processes  such  as  extracellular  matrix  reorganization 
and  inhibition  of  apoptosis  suggests  that  cell  cycle-dependent  transcriptional  control 
influences  the  execution  of  many  biological  programs  not  obviously  related  to  mitosis.  In 
marked  contrast  to  the  limited  changes  observed  following  deletion  of  specific  transcription 
factors26,48;  the  regulatory  fluctuations  observed  here  and  in  profiling  of  other  general 
cellular  processes  involve  broad  functional  classes  of  genes  and  multiple  members  of  gene 
families*  2,49.  Analogous  to  quantitative  genetic  phenomena,  certain  biological  activities 
may  depend  on  the  aggregation  of  transcriptional  effects  whose  significance  is  difficult  to 
measure  on  an  individual  basis. 

The  functional  classifications  applied  in  this  study  represent  only  a  first  attempt  to 
formally  structure  public  domain  knowledge  about  human  genes  for  the  purpose  of 
systematically  analyzing  experimental  genomic  data.  Generation  of  increasingly  refined 
representations  of  such  information  should  enable  more  a  far  more  sophisticated 
understanding  of  this  and  other  transcriptional  data  sets.  Certainly,  many  biological 
activities  may  be  mediated  through  a  limited  number  of  transcriptional  targets  and  will 
therefore  not  meet  the  statistical  thresholds  inherent  to  this  approach.  However,  given  the 
considerable  size  and  proliferation  of  expression  data  sets,  functional  evaluations  now 
represent  the  rate-limiting  step  to  deriving  meaning  from  these  studies.  The  methods 
demonstrated  here  provide  a  first  solution  to  this  growing  problem.  Because  of  the  scope  of 
this  data  set,  more  conclusions  remain  to  be  drawn  than  can  be  described  in  this  manuscript. 
We  encourage  the  reader  to  directly  access  the  data  via  the  Internet  or  FTP. 

As  for  the  yeast  DNA  damage  inducible  genes,  we  have  found  a  large 
number  of  genes  that  were  inducible  in  a  Mecl  and  Dunl-dependnet  manner.  Several  sets 
of  these  inducible  genes  are  represented  by  general  oxidative  stress  or  by  protein  unfolding. 
We  observed  a  large  number  of  proteins  involved  in  protein  folding  as  well.  These  will  serve 
as  a  point  of  comparison  for  the  mammalian  genes  isolated  now  and  in  the  future  to 
determine  the  similarities  between  organismal  responses  to  DNA  damage. 

Methods 

Synchronization  of  Cells 

Human  cell  lines  used  in  this  study  were  derived  from  human  foreskin  samples^. 
Cells  were  grown  in  Dulbecco’s  modified  Eagle’s  medium  supplemented  with  10%  fetal 
bovine  serum  and  antibiotics  (all  from  GIBCO).  Cells  were  grown  in  150  cm2  flasks  to  25- 
30%  confluence  and  synchronized  using  a  double  thymidine  block  protocol^,  with  a  first 
block  for  16  h,  a  10  h  release,  and  a  second  block  for  15  h.  The  final  concentration  of 
thymidine  used  in  the  block  medium  was  2  mM.  Approximately  2  x  107  cells  were 
harvested  at  each  time  point.  About  10%  of  cells  from  each  time  point  were  fixed  in  ethanol 
and  stained  with  propidium  iodide  for  flow  cytometry  analysis,  and  the  remainder  were 
immediately  frozen  at  -80  °C  for  subsequent  isolation  of  RNA. 
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UV  and  MMS  treatment 

Cells  were  grown  asynchronously  and  passaged  before  reaching  80%  confluency.  At 
the  last  passage,  cells  were  combined  and  split  equally  to  18  150  cm2  dishes  and  grown 
overnight  (~16  h).  Media  was  then  removed  from  each  dish  and  stored  separately.  All 
dishes  were  then  rinsed  once  with  IX  PBS.  Six  dishes  were  treated  with  UV  irradiation  (40 
J/m2)  and  six  dishes  were  treated  with  MMS  at  a  final  concentration  of  0.015%  (v/v).  The 
remaining  six  dishes  were  not  subjected  to  any  treatment.  Media  was  then  added  back  to  all 
dishes  and  all  samples  were  incubated  another  4  h  before  harvest.  For  harvest,  cells  were 
trypsinized  for  removal  from  plates,  immediately  flash  frozen  in  liquid  nitrogen,  and  stored 
at  -80  °C . 

Isolation  and  Hybridization  of  mRNA 

Total  RNA  was  isolated  and  biotin-labeled  as  previously  described  10.  Samples  were 
hybridized  to  a  total  of  five  arrays  containing  approximately  1.4  x  106  oligonucleotide 
probes  to  approximately  40,000  human  genes  and  ESTs  (Affymetrix,  Santa  Clara). 
Hybridizations  were  carried  out  at  42  °C  for  14  to  16  h  with  mixing  on  a  rotisserie  at  60 
rpm.  Following  hybridization,  the  solutions  were  removed,  the  arrays  were  rinsed  with  6X 
SSPE-T  (0.9  M  NaCl,  60  mM  NaH2P04,  6  mM  EDTA,  0.005%  Triton  X-100  adjusted  to 
pH  7.6),  rinsed  with  0.5X  SSPE-T  (75  mM  NaCl,  5  mM  NaH2P04,  0.5  mM  EDTA, 
0.005%  Triton  X-100  adjusted  to  pH  7.6),  and  incubated  with  0.5X  SSPE-T  at  42  °C  for  15 
minutes.  Following  washing,  hybridized  biotinylated  cRNA  was  fluorescently  labeled  by 
incubation  with  2  pg/ml  streptavidin-phycoerythrin  (Molecular  Probes,  Eugene,  OR)  and  1 
mg/ml  acetylated  BSA  (Sigma,  St.  Louis,  MO)  in  6X  SSPE-T  at  42  °C  for  10  min. 
Unbound  streptavidin-phycoerythrin  was  removed  by  rinsing  at  room  temperature  prior  to 
scanning.  The  arrays  were  read  at  a  resolution  of  3  pm  using  a  specially  designed  confocal 
scanner  (Affymetrix,  Santa  Clara,  CA)  as  described  previously^. 

Expression  Data  Analysis 

Data  were  normalized  between  samples  by  setting  the  mean  hybridization  signal  for 
each  sample  equal  to  1000  arbitrary  units  and  the  standard  deviation  to  1500  arbitrary  units 
after  discarding  the  top  2%  and  the  bottom  2%  of  all  signals.  The  hybridization  signal  for 
each  gene  or  EST  was  then  normalized  to  a  unit  standard  deviation  with  a  mean  of  zero. 
Genes  having  similar  intensity  patterns  across  multiple  time  points  were  detected  by 
clustering  using  the  Pearson  correlation  coefficient  as  the  metric  of  similarity.  Genes  were 
placed  into  specific  clusters  by  identifying  relevant  expression  patterns  and  matching  the 
expression  profile  of  each  gene  to  a  pattern  using  a  correlation  coefficient  threshold. 

The  abundance  of  cell  cycle  phase-specific  transcripts  is  expected  to  fluctuate 
periodically.  We  identified  periodic  patterns  using  supervised  and  unsupervised  clustering 
methods.  Supervised  clustering  was  based  on  comparisons  to  series  of  sine  waves  offset  bv 
consecutive  time  intervals,  similar  to  previous  analyses  of  cell  cycle  expression  data®. 
Typically,  we  demanded  that  a  gene  be  correlated  to  one  of  the  sine  wave  patterns  by  a 
Pearson’s  correlation  coefficient  of  at  least  0.70  over  24  samples.  However,  many  known 
cell  cycle-regulated  transcripts  displayed  patterns  that  did  not  conform  to  simple  sine  waves. 
Therefore,  we  also  searched  for  common  patterns  in  the  data  by  conducting  unsupervised 
clustering  between  all  genes,  searching  for  groups  of  at  least  four  transcripts  that  exhibited  a 
correlation  coefficient  of  at  least  0.90.  Twenty-two  such  groups  were  identified,  of  which 
five  were  found  to  be  cell  cycle-regulated.  These  groups  were  then  averaged  and  used  as 
seed  patterns  for  clustering  against  all  genes,  using  a  correlation  coefficient  of  at  least  0.70. 
Groupings  from  supervised  and  unsupervised  clustering  with  identical  peak  times  were 
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merged,  resulting  in  the  clusters  described  in  the  text.  More  than  90%  of  these  transcripts 
were  scored  as  detectable  in  more  than  a  quarter  of  the  timepoints. 

A  database  of  functional  information  for  more  than  2400  human  genes  was 
generated  by  manual  placement  of  genes  into  functional  categories  using  information  from 
Medline  and  other  publicly  available  information  sources.  The  database  contains  160 
controlled-vocabulary  categories  that  are  classified  hierarchically.  The  database  schema  is 
similar  to  that  of  the  MIPS  yeast  database  (www.mips.org),  but  also  contains  functions 
relevant  to  multicellular  organisms.  This  database  is  commercially  available  through 
Molecular  Applications  Group  (www.mag.com). 

The  binomial  distribution  function  was  applied  using  the  Stingray  Expression 
Analysis  software  package  (Molecular  Applications  Group,  Palo  Alto,  CA).  Briefly,  if  the 
number  of  genes  in  an  expression  cluster  with  a  given  biological  function  was  less  than  or 
equal  to  the  number  expected  by  random  chance,  then  the  P  value  was  calculated  using  the 
formula: 


n 

x 


px(i-p) 


n-x 


where  n  is  the  number  of  genes  in  a  given  cluster,  x  is  the  observed  number  of  genes  found 
in  both  a  given  cluster  and  a  given  functional  category,  and  p  =  the  overall  frequency  at 
which  genes  of  a  given  functional  category  are  found  in  the  genes  represented  on  arrays.  If 
the  number  of  genes  in  an  expression  cluster  with  a  given  biological  function  was  greater 
than  the  number  expected  by  random  chance,  then  the  P  value  was  calculated  using  the 
formula: 


Table  1  Human  genes  expressed  maximally  in  G2/M. 


Accession  Number  Gene  function 


X66945_at  3 

X54942_at  3 

X79201_at  3 

Z36714_at  3 

X87613_at  3 

Z49989_at  3 

L25876_at  3 

M62994_at  3 

U47635_at  3 

X14850_at  3 

U14518_at  3 

M25753_at  3 

U73379_at  3 

M96803_at  3 


H.sapiens  N-sam  mRNA  for  fibroblast  growth  factor  receptor 

H.sapiens  ckshs2  mRNA  for  Cksl  protein  homologue 

H.sapiens  mRNA  for  SYT 

H.sapiens  mRNA  for  cyclin  F 

H.sapiens  mRNA  for  skeletal  muscle  abundant  protein 

H.sapiens  mRNA  for  smoothelin 

Homo  sapiens  protein  tyrosine  phosphatase  (CIP2)mRNA, 

Homo  sapiens  thyroid  autoantigen  (truncated  actin-binding  protein) 
Human  D13S824E  locus  mRNA,  complete  cds 
Human  H2A.X  mRNA  encoding  histone  H2A.X 
Human  centromere  protein-A  (CENP-A)  mRNA,  complete  cds 
Human  cyclin  B  mRNA,  3'  end 

Human  cyclin-selective  ubiquitin  carrier  protein  mRNA, 

Human  general  beta-spectrin  (SPTBN1)  mRNA,  complete  cds 
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Accession  Number 
M86699_at  3 

U37426_at  3 

D26361_at  3 

D31885_at  3 

D43948_at  3 

X51688_at  3 

D78514_at  3 

U63743_at  3 

U28386_at  3 

U05340_at  3 

U01038_at  3 

M90657_at  3 

M91670_at  3 


Gene  function _ 

Human  kinase  (TTK)  mRNA,  complete  cds 

Human  kinesin-like  spindle  protein  HKSP  (HKSP)  mRNA, 

Human  mRNA  for  KIAA0042  gene,  complete  cds 
Human  mRNA  for  KIAA0069  gene,  partial  cds 
Human  mRNA  for  KIAA0097  gene,  complete  cds 
Human  mRNA  for  cyclin  A 

Human  mRNA  for  ubiquitin-conjugating  enzyme,  complete  cds 

Human  mitotic  centromere-associated  kinesin  mRNA,  complete  cds. 

Human  nuclear  localization  sequence  receptor  hSRPl  alpha  mRNA, 

Human  p55CDC  mRNA,  APC  activator 

Human  pLK  mRNA,  complete  cds 

Human  tumor  antigen  (L6)  mRNA,  complete  cds 

Human  ubiquitin  carrier  protein  (E2-EPF)  mRNA,  complete  cds 


Table  2.  Human  genes  maximally  expressed  during  S  phase. 

Accession  Number  Gene  function 


X68277_at  7 

X66087_at  7 

X89750_at  7 

Z46967_at  7 

X62535_at  7 

Z24725_at  7 

M77140_at  7 

Z24727_at  7 

U62015_at  7 

Z46629_at  7 

J02854_at  7 

U73960_at  7 

U44975_at  7 

U27655_at  7 

U27768_at  7 

U59752_at  7 

U28049_at  7 

X16416_at  7 

M92934_at  7 

U04636_at  7 

U18300_at  7 

U81607_at  7 

X52142_at  7 

D90209_at  7 

X52599_at  7 

X58377_at  7 

L08246_at  7 

U08021_at  7 

U40490_at  7 

M12174_at  7 

M62831_at  7 

J03764_at  7 


H.sapiens  CL  100  mRNA  for  protein  tyrosine  phosphatase 
H.sapiens  a-myb  mRNA 
H.sapiens  mRNA  for  TGIF  protein 
H.sapiens  mRNA  for  calicin  (partial) 

H.sapiens  mRNA  for  diacylglycerol  kinase 

H.sapiens  mitogen  inducible  gene  mig-2,  complete  CDS 

H.sapiens  pro-galanin  mRNA,  3'  end 

H.sapiens  tropomyosin  isoform  mRNA,  complete  CDS 

Homo  sapiens  Cyr61  mRNA,  complete  cds 

Homo  sapiens  SOX9  mRNA 

Human  20-kDa  myosin  light  chain  (MLC-2)  mRNA,  complete  cds 

Human  ADP-ribosylation  factor-like  protein  4  mRNA,  complete  cds 

Human  DNA-binding  protein  CPBP  (CPBP)  mRNA,  partial  cds 

Human  RGP3  mRNA,  complete  cds 

Human  RGP4  mRNA,  complete  cds 

Human  Sec7p-like  protein  mRNA,  partial  cds 

Human  TBX2  (TXB2)  mRNA,  complete  cds 

Human  c-abl  mRNA  encoding  pi 50  protein 

Human  connective  tissue  growth  factor,  complete  cds 

Human  cyclooxygenase-2  (hCox-2)  gene,  complete  cds 

Human  damage-specific  DNA  binding  protein  p48  subunit  (DDB2) 

Human  gravin  mRNA,  complete  cds 

Human  mRNA  for  CTP  synthetase  (EC  6. 3.4.2) 

Human  mRNA  for  DNA  binding  protein  TAXREB67 
Human  mRNA  for  beta  nerve  growth  factor 
Human  mRNa  for  adipogenesis  inhibitory  factor 
Human  myeloid  cell  differentiation  protein  (MCL1)  mRNA 
Human  nicotinamide  N-methyltransferase  (NNMT)  mRNA, 

Human  nicotinamide  nucleotide  transhydrogenase  mRNA, 

Human  ras-related  rho  mRNA  (clone  6),  partial  cds 
Human  transcription  factor  ETR101  mRNA,  complete  cds 
Human,  plasminogen  activator  inhibitor- 1  gene,  exons  2  to  9 
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(6)  Key  Research  Accomplishments 

A)  Identification  of  all  known  DNA  damage  inducible  genes  in  the  yeast  S.  cerevisiae. 
Described  in  previous  report 

B)  Determination  of  cell  cycle  regulated  genes  in  human  diploid  fibroblasts,  40,000 
transcripts  analyzed. 

C)  Determination  of  DNA  damage-inducible  genes  in  human  diploid  fibroblasts,  6,800 
transcripts  analyzed. 

(7)  Reportable  Outcomes. 

Paper  submitted  to  Molecular  Biology  of  the  Cell,  in  review.  A  second  paper 
published  in  Nature  Genetics 

Cho,  R.J.,  Huang,  M.,  Campbell,  M.J.,  Dong,  H.,  Steinmetz,  L.,  Sapinoso,  L.,  Hampton,  G., 
Elledge,  S.J.,  Davis,  R.W  and  Lockhart,  D.J.  (2001)  Transcriptional  regulation  and  function 
during  the  human  cell  cycle.  Nat.  Genet.  27:  48-54. 

(8)  Conclusions 

So  far  our  work  has  established  that  a  large  number  of  human  genes  are  controlled 
in  a  cell  cycle  regulated  manner.  Many  processes  previously  not  known  to  be  cell  cycle 
regulated  were  found  to  be  regulated  in  this  manner.  Furthermore,  many  cell  cycle  regulated 
genes  are  also  inducible  by  DNA  damage. 
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